HFOSR  -TR-.  76  - I '02  7 


. Final '"Report 
on  the 


30  •< 

rh ; 
io 


«5l 

jO 

■ iso 


REGIONALIZATION  OF  THE  ARCTIC  REGION,  SIBERIA 
AND  EURASIAN  CONTINENTAL  AREA 


Sponsored  by 


Advanced  Research  Projects  Agency 
ARPA  Order  No.  1827 


D D C 

SEP  13  1976 


Principal  Investigator: 


Program  Manager: 


Contract  Number: 


Amount  of  Contract:  / 

Effective  Dates  of  Contract 
Program  Code: 


Professor  Leon  Nnopoff 
(Telephone  No:  213-825-1885) 

Mr.  Wil I iam  J.  Best 
(Telephone  No:  202-694-5456) 

Air  Force  Contract  F44620- 
73-C-0048 


$164,986 

March  1,  1973  to  May  31,  1976 
4 FI  0 


i :i" . , • - ; 

Institute  of  Geophysics  and 
University  of  California 

r StlpSSfc  «— i JjZZB 

I , i? 13  ' ■ 

■ D.  bloss  W0-i3  (7bJ_ 

«« 0«ioer 


Planetary  Physics 
, Los  Angeles 


>■ 


table  of  contents 


Technical  Report  Summary  . 


Technical  Report 

I.  Purposes  of  project  

"•  Review  of  scientific  background 



HI.  Objectives,  methods  and  results 
IV.  Bibliography 
V.  Figures  . . 


utasm  ror “ 

****  WJilfe  Sstiion  [y 

. MI  Section  □ 

WAHHOIWCEB  a 

JI/STIflCAHON , 


(**»«! 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  THE  BEST 
QUALITY  AVAILABLE. 

COPY  FURNISHED  CONTAINED 
A SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


\ 


TECHNICAL  REPORT  SUMMARY 


In  this  investigation,  our  first  purpose  is  the  region- 
alization of  the  Arctic  region,  Siberia  and  the  Eurasian  con- 
tinental area  using  seismic  surface  waves.  This  regionaliza- 
tion will  determine  the  structural  properties  of  the  upper  few 
hundred  kilometers  of  the  earth,  and  the  variation  of  these 
properties  from  one  subregion  to  another  in  the  area  under  in- 
vestigation. Once  the  structural  properties  in  the  various 
regions  have  been  obtained,  our  second  purpose  is  to  apply  op- 
timized computer  techniques  to  the  computation  of  accurate 
theoretical  seismograms  for  any  hypothetical  type  of  source 
located  anywhere  within  this  continental  area.  These  theoret- 
ical seismograms  can  be  applied  directly  to  the  discrimination 

problem  by  comparing  sei smograms, computed  for  both  earthquakes 
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Regional ization.  The  first  technical  problem  to  be  dealt 
with.,  which  is  complete,  is  the  computation  of  experimental  phase 
velocity  curves  for  surface  waves  traversing  all  of  the  regions 
under  investigation. 

The  phase  velocity  curves  were  obtained  with  the  single- 
station phase  velocity  method,  which  is  described  in  the  main 
text  of  this  report.  Briefly,  the  method  involves  scanning  our 
microfilm  library  of  World  Wide  Standardized  Sei smographic  Net- 
work (WWSSN)  records  for  earthquakes  which:  (1)  occurred  within, 
or  on  the  perimeter  of,  the  area  of  interest,  (2)  produced  good 
long-period  surface  wave  records  at  WWSSN  stations  for  which  the 
epicenter-to-station  lines  lie  within  the  regions  being  investi- 
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gated,  and  (3)  generated  good,  large  recordings  at  a sufficient 
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number  of  stations  to  ensure  an  accurate  fault  plane  solution. 

When  a suitable  earthquake  has  been  found,  an  extensive  data  pro- 
cessing and  data  reduction  system  is  applied  to  transform  the  data 
into  phase  velocity  curves  for  each  of  the  selected  epicenter- 
to-station  lines  through  the  area  of  interest. 

The  complete  set  of  phase  velocity  curves  (for  fundamental- 
mode  Rayleigh  waves  recorded  by  the  WWSSN  Instruments),  which  will 
be  required  for  the  regionalization  of  the  Arctic  region,  Siberia, 
end  the  Eurasian  continental  area,  consists  of  about  50  dispersion 
curves.  Based  on  this  data  set,  we  have  just  obtained  our  first 
really  successful  regionalization  of  the  structures  in  the  area 
under  Investigation.  The  computer  programs,  upon  which  the  inver- 
sion of  the  data  for  structural  information  is  based,  were  first 
used  successfully  in  the  regionalization  of  the  Pacific  Ocean  area 

IKausel,  19/2;  Leeds,  1973;  Kausel,  Leeds  and  Knopoff,  1974;  Leeds, 
Knopoff  and  Kausel,  1974), 

Our  attack  on  the  Inversion  problem  has  been  a three-part 
effort:  (1)  the  regionalization  of  Siberia  and  the  Eurasian  conti- 
nental area,  using  only  the  data  for  paths  limited  to  these  regions; 
-<2>  the  regionalization  of  the  Arctic  area,  again,  using  only  the 
data  for  paths  limited  to  this  region;  and  (3)  the  final  regionali- 
zation of  Siberia,  the  Eurasian  continental  area,  and  also  the 
Arctic  region,  using  the  entire  set  of  data  we  have  available,  in- 
cluding many  profiles  having  both  continental  and  oceanic  segments 
between  epicenter  and  station.  We  have  completed  the  first  success- 
ful inversion  process  relative  to  part  (1),  and  reduced  the  data 
from  part  (2)  to  obtain  the  phase  velocity  for  the  subregion  of 
major  interest  in  this  area;  the  execution  of  part  (3)  is  planned 


3. 


ar%5*sSirt(i*W»^  • - — •*■-*«  T/.w:  • u"cx>uriy&*'-s-nr--‘-f>**'  /’-•  •■' 


for  our  next  contract  year. 

From  our  preliminary  regionalization  work  on  the  conti- 
nental areas,  we  can  draw  a few  conclusions  already.  First, 
we  note  that  the  properties  of  much  of  the  Asiatic  upper  man- 
tle, averaged  over  long  trans-Eurasian  distances,  are  highly 
consistent  with  those  of  young  stable  regions  observed  else- 
where (Knopoff,  1972;  Fouda,  1973;  Knopoff  and  Fouda,  1974). 
Second,  the  Tibetan  plateau  has  an  extremely  thick  crust, 

, perk  s as  great  as  80  km  from  surface  to  Moho.  Third, 

the  A.pide  folded  belt  of  Iran  and  Turkey  has  an  extraordina- 
rily wel 1 -devel oped  low-velocity  channel  in  the  mantle,  of  re- 
markable contrast  to  the  lid  above.  The  principal  problem  in 
the  inversion  has  been  the  construction  of  the  boundaries  to 
the  geological  provinces,  for  which  only  incomplete  information 
is  found  in  the  literature.  This  difficulty  is  what  has  led  us 


to  the  idea  of  repeating  parts  (1)  and  (2)  of  our  scheme  of  re- 
gionalization. The  detailed  (preliminary)  results  of  our  inver- 

slon  are  given  on  pages  21-38. 

Theoretical  Seismograms.  The  second  phase  of  our  con- 
tract work  ~ computing  realistic  theoretical  seismograms  at 
the  WWSSN  stations,  and  especially  at  the  new  high-gain,  wide- 
band installations  (Pomeroy  eJLii*  * 1 969;  Molnar  » 1 969, 

Savlno,  McCarry  and  Hade. 1972)  at  Kongserg,  Norway  (ICON.)  * Toledo, 
Spain  (TOL),  Eilat,  Israel  (EIL)»  Chiengmai , Thailand  ( CHG ) , and 
Matsushiro,  Japan  (MAT),  for  arbitrary  seismic  and  explosive 
sources  within  the  area  under  invest! agation  - is,  of  course, 
dependent  upon  the  availability  of  the  final  regionalization 
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from  the  first  section  of  our  investigation.  To  prepare  for 
this  work,  we  have  concentrated  upon  improving  our  existing 
algorithms  and  computer  programs  for  generating  theoretical 
seismograms  for  highly  realistic  models  of  the  earth.  In 
this,  we  have  succeeded  to  the  extent  that  only  the  problem 
of  the  lateral  heterogeneity,  over  the  area  of  investigation, 
remains  as  an  important  consideration  in  the  further  develop- 
ment of  our  system  for  computing  these  theoretical  time  series. 
These  statements  apply  to  the  LoveJ  or  torsional  waves  excited 
by  earthquakes  and  explosions;  the  parallel  development  of  our 
algorithms  and  computer  programs  for  Rayleigh,  or  spheroidal 
waves  on  a spherical,  gravitating  earth  is  currently  one  of  the 
main  efforts  under  this  contract. 

The  work  published  out  of  our  section  has  become  quite 
voluminous  (see  following  Technical  Report)  in  the  subjects  -- 
dispersion  computations,  attenuation  computations,  structural 
transformations,  eigenfunction  characteristics,  effects  of 
sphericity,  point-source  response  computations  — upon  which 
efficient  construction  of  accurate,  multimode  theoretical  seis- 
mograms for  realistic  models  of  the  earth's  structure  are  based. 

Our  original  work  with  theoretical  seismograms  concerned 
the  multimode  surface  wave  phases  Lg  (Knopoff,  Schwab  and  Kau- 
sel,  1973;  Knopoff  et  al_. , 1974)  and  Sa  (Schwab  e_t  al_.  , 1974). 
Relative  to  discrimination  studies,  however,  it  is  desirable  to 
have  body  wave  phases’  also  available  on  the  theoretical  seismo- 
grams; this  is  one  direction  of  our  current  work;  the  other  is 
improvement  in  efficiency,  control  of>accuracy,  and  increased 


power  and  flexibility  of  the  computational  algoritms  and  com- 
puter programs.  Our  first  successes  in  generating  complete 
seismograms,  l.e.,  body  and  surface  waves  on  the  same  record, 
were  reported  by  Nakanlshi  et  al . , (1  975),  Nakanishi,  Schwab 
and  Kausel  (1976),  and  Nakanlshi,  Schwab  and  Knopoff  (1976). 

This  work  includes  the  summation  of  eight  modes,  and  the  com- 
putations are  applied  to  a continental  structure;  our  first 
multimode  time  series  for  an  oceanic  structure  (Kausel,  Schwab 
and  Mantovani,  1976)  also  forms  part  of  this  research.  Multi- 
mode  theoretical  seismograms  containing  body  wave  phases  have 
been  generated  in  the  past  (Sato,  Usami  and  Landisman,  1968); 
however,  these  time  series  were  limited  to  ultralong  periods. 

The  Important  point  of  our  present  results  is  that,  owing  to  the 
new  high  efficiency  of  our  algorithms,  we  have  been  successful 
in  extending  the  period  content  of  the  theoretical  seismograms 
down  through  the  range  covered  by  the  long-period  instruments 
of  the  WWSSN.  Thus,  we  can  generate  the  theoretical  time  series 
which,  for  the  first  time  we  believe,  permits  us  to  begin  think- 
ing about  the  comparison  of  theoretical  seismograms  with  the 
entire  records  obtained  at  the  WWSSN  and  the  new  high-gain,  wide- 
band installations.  Our  current  level  with  these  theoretical- 
numerical  techniques  permits  the  generation  of  time  series  with 
a combination  of  eleven  modes,  and  for  the  arbitrary  specifica- 
tion of  source  depth  as  well  as  the  other  focal  parameters.  The 
results  of  this  level  of  our  system  are  given  by  Mantovani  ert  al_.  (1976a, 
1976b).  Our  immediate  goal  in  this  work  is  the  extension  to 
the  combination  of  fifteen,  and  finally  twenty,  modes  to  form 
these  theoretical  time  series. 
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technical  report 

1 Purposes  of  this  project 

H * n^niprt  is  to  determine  the 

The  first  purpose  of  our  prooect 

t.  , ,nd  the  variation  of  these  properties, 
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to  the  long-period  surface-wave  1nformation 
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.structural  properties  of  Siberia  r 

51 Deri a , the  Eurasian  continental 

area,  and  the  Arctic  region. 

Review  of  scientific  background 

Mionalization.  Since  the  single-station  surface 
wave  method  allows  al,  stations  to  be  located  at  the  edge  of 
the  region  under  investigation,  this  technique  is  ideal  for 
the  proposed  study.  Brune  et  al..  (I960)  were  the  first  to 
describe  the  single-station  method.  Knopoff  and  Schwab  (1968) 
corrected  and  extended  this  description  to  take  into  account 
the  frequency  dependence  of  the  apparent  initial  phase  of  the 
source.  Based  on  the  Thomson  0 95 0 ) -Has ke 1 1 (,953)  matrix  for- 
mulation for  surface  wave  propagation,  Harkrider  (1964)  devel- 
oped the  theory  for  treating  the  surface  wave  response  to  buried 
point-source  singlets.  Ben-Menahem  and  Toksbz  (1963)  then  de- 
veloped the  formalism,  in  terms  of  the  singlet  response,  for 
representing  the  displacement  field  of  an  arbitrary  force  sys- 
tem in  a multilayered  medium;  and  Ben-Menahem  and  Harkrider 
(1964)  applied  this  formalism  to  obtain  the  surface  wave  re- 
sponse to  point-source  couples  and  double  couples.  Harkrider 
(1970)  later  made  certain  corrections  and  Improvements  in  these 
last  results.  The  far-field  response  to  real.istic,  displacement- 
dislocation  faulting  was  shown  to  be  equivalent  to  the  point- 
source,  double-couple  force  replacement  in  an  unfaulted  medium 
(Maruyama,  ,963;  Burridge  and  Knopoff,  ,964),  and  the  single- 
station method  was  on  firm  theoretical  ground  for  flat,  multi- 
layered media.  By  means  of  transformation  techniques  (Biswas 
and  Knopoff,  1970;  Schwab  and  Knopoff,  1972),  it  is  possible  to 
convert  point-source  programs  for  a flat  structure  into  those 
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useful  for  treating  $r  ces  A cnu  . 
a . c u 1n  a sPhencal  structure  (Kausel 

C wab,  1 973).  Ti  . programs  for  the  surface,  mve 

sphere  '■  ■ heterogeneous 

sphere  are  now  a part  of  our  program  1 1 bpa ry . We  are  „s  . 

entire "e9c’tand  °Ptl"”Zln9  Pr09r“S  as  P»«  •'  the  pres- 

-«sper,I.TdB.1rwrTd^\n,dnVtr#CtUr*1  1nf°r"ti0n  fr°n' 

nr  . sc“ssed  by  Knopoff  (1  961.1962),  and  the 

Prmn9  °f  thS  procedure  appears  to  have 

been  carried  out  by  Dorman  and  Ewing  (,962).  These  papers  de- 

;;;ed  auempts  to  ^ f r 
r~’  — — were  mainly  concerned 

within  th  ^ 6 ^ StrUCtUral  "de,S  “hiCh  fU  thC  data  *o  - 

he  experimental  accuracy  ( Kei, is-Borok  and  Yanovskaya, 

i h 7"  ,968,'969)-  ™S  "*  «•  P~ming,  un- 

direction  of  Knopoff  and  KeiTis-Borok,  of  the  Hed9eho9 
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isfy  ;h  TS  te  'qUe  inV0lV6S  Seek,'n9  StrUCtU-’  »t- 

• . ata  t0  WUInn  the  exPerlmental  accuracy,  by  means  of 

-twise  search  throughout  a multidimensional  parameter  space 
, package  is  efficient  enough  for  routine  application  to  sur- 
face wave  dispersion  data.  «e  have  now  introduced  variational 
Parameters  into  the  procedure,  which  has  resulted  in  a signifi- 
cant  increase  in  computational  efficiency.  A description  of 

and  e"S-1PVOlVed  1n  thlS  lrnP™vement  is  given  by  Wiggins  (1  972) 

"d  dackson  (.972),  the  technigue  takes  into  account  the  redun- 
dancy in  the  data  as  well  as  numerical  instabilities  in  the 
solution. 


The  single-station,  surface  wave  regionalization  of  an 
area  is  based  on  the  assumption  of  great-circle  propagation 
paths,  and  the  assumption  that  the  phase  travel  time  from  epi- 
center to  station  is  just  the  sum  of  the  travel  times  through 
the  subdivisions  of  the  laterally  heterogeneous  region  (Knopoff, 

1969) .  This  permits  the  regionalization  to  be  expressed  at 
each  frequency  as  a system  of  inhomogeneous  equations;  for  each 
path,  the  total  travel  time  is  the  inhomogeneous  term,  the  dis- 
tances through  the  subdivisions  are  the  coefficients,  and  the 
slownesses  in  the  subdivisions  are  the  unknowns.  Santo  [equa- 
tion (1),  1961b]  applied  this  technique  to  the  Pacific,  using 
group  travel  times  and  group  slownesses  but,  due  to  considera- 
tions involving  uniqueness  of  inversion  (Pilant  and  Knopoff, 

1970) ,  we  prefer  phase  travel  time  and  phase  slownesses.  The 
solution  of  this  system  of  equations  yields  the  experimental 
phase  slownesses  associated  with  each  of  the  subdivisions,  in 
principle;  however,  the  above  procedure  represents  a version 
of  the  computer  program  which  is  not  pedagogi cal ly  acceptable. 
Instead,  we  prefer  to  consider  the  model  parameters  as  the  pri- 
mary unknowns  in  the  inversion  and  to  use  the  phase  velocities 
therefrom. 

The  initial,  large  scale,  lateral  regionalization  of  a 
portion  of  the  earth  by  means  of  single-station  surface  wave 
studies  was  performed  by  Santo  (1960a;  1960b;  1961a;  1961b)  in 
an  investigation  of  the  Pacific  Ocean  basin.  The  bandwidth  of 
his  recordings  and  his  data  processing  techniques  limited  al- 
most all  of  his  dispersion  results  to  periods  less  than  about 
40  seconds;  Santo's  regionalization  did  not  yield  structural 
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information  below  the  crust  and  lid  of  the  mantle.  Our  current 
data  processing  techniques  (Biswas,  1971),  and  the  availability 
of  a nearly  complete  library  of  long-period  seismograms  from 
the  WWSSN  stations,  have  made  it  possible  to  extend  the  lateral 
regionalization  to  a depth  of  about  250  km.  The  first  applica- 
tion of  our  single-station  regionalization  techniques  was  an 
investigation  of  the  upper  mantle  structure  in  the  Pacific 
Ocean  Basin  (Leeds,  1973). 

To  exhibit  the  potential  power  of  this  technique,  a 
brief  summary  of  our  work  with  the  structure  in  the  Pacific 
Ocean  Basin  will  be  useful:  The  results  of  the  work  of  Kausel 

(1972)  showed  an  apparently  continuous  gradation  of  Rayleigh- 
wave  phase  velocities  from  low  values  observed  on  paths  close 
to  the  East  Pacific  Rise  to  relatively  high  values  on  paths  in 
the  oldest  parts  of  the  lithosphere.  Leeds  (1973)  increased 
the  body  of  data  and  performed  the  first  successful  inversion 
of  the  data,  assuming  that  regional  variations  in  structure 
were  correlated  with  lithospheric  age;  that  study  took  into 
account  bathymetry  and  available  crustal  information.  The  num- 
ber of  degrees  of  freedom  in  the  data  was  found  to  be  remark- 
ably small,  despite  the  large  number  of  phase-velocity  determi- 
nations. The  inversion  led  to  the  conclusion  that  the  data  set 
did  not  permit  one  to  obtain  detailed  information  concerning 
the  bottom  of  the  low-velocity  channel.  Furthermore,  the  lith- 
osphere increased  in  thickness  monotonically  with  age: 
at  the  ridge  crest,  the  lithosphere  has  almost  zero  thick- 
ness. This  is  thus  the  explanation  for  the  observation 
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that  the  phase  velocities  change  systematically  with  distance 
from  the  ridge  crest.'  The  model  of  lithospheric  thickness  as 
a function  of  spreading  age  is  consistent  with  the  model 
wherein  the  1 i thosphere-a sthenosphere  interface  is  at  the 
solidus  for  wet  peridotite.  The  results  of  the  Kausel  and 
the  Leeds  Ph.D.  theses,  as  well  as  some  later  conclusions, 
now  appear  in  the  literature:  Kausel,  Leeds  and  Knopoff  (1974), 
Leeds,  Knopoff  and  Kausel  (1974),  and  Leeds  (1975).  These  single- 
station regionalization  techniques  are  being  refined  for  use  in 
the  present  project. 

Upper  mantle  studies  in  Eurasia  appear  to  have  suffered 
due  to  the  Soviet  seismologists'  lack  of  efficient  long-period 
instruments.  Also,  their  extensive  program  in  deep  seismic 
sounding  has  undoubtedly  served  to  focus  interest  upon  crustal 
•studies,  to  the  detriment  of  sub-Moho  investigations.  The- 
review  by  Arkhangel 1 skaya  (1960)  discusses  the  foreign  and  do- 
mestic surface  wave  studies  performed  in  Eurasia  prior  to  1960. 
The  early  Soviet  studies  discussed  in  this  review,  as  in  the 
later  studies  to  be  mentioned  here,  were  limited  to  short-period 
investigations  --  usually  less  than  40  seconds  --  and  were  con- 
cerned mainly  with  determining  crustal  properties.  Of  the  early 
Soviet  studies,  which  more  or  less  paralleled  the  work  in  the 
west,  the  investigations  of  Savarensky  and  Ragimov  (1958;  1959), 
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Savarensky,  Solov'eva  and  Shechkov  (1959),  and  Savarensky  and 

S<kharU,idZe  °953)  demonstrate  this  predominant  inter- 

est in  crustal  features.  The  subsequent  work  by  Popov  (1960), 

Shechkov  (1961),  Savarensky  and  Shechkov  (1961),  Shechkov  and’ 
Solov'eva  (1961),  and  Shechkov  (1964)  indicated  that  the  Soviet 
surface  wave  investigations  would  remain  focused  on  crustal 
studies,  and  ideed,  this  has  proven  to  be  the  case  up  through 
the  most  recent  Soviet  surface  wave  literature  available  to  us 
(Savarensky  and  Peshkov,  1968;  Sikharulidze  and  Makharadze,  1968; 
Savarensky  et  aJL  , 1 969  , Shechov,  1 970). 


All  of  these  crustal  studies,  as  well  as  Santo's  (1962, 
1965a;  1965b)  short-period  regionalization  efforts  should  be  val- 
uable in  assisting  the  determination  of  which  shallow  structures 
to  combine  with  our  trial  parameterizations  for  the  deep  struc- 
tures in  the  inversion  portion  of  the  present  investigation. 

Ih^«M«l^eismocm^.  Relative  to  the  discrimination 
problem,  probably  the  most  important  feature  in  the  calculation 
of  theoretical  seismograms  which  requires  improvement  over  pre- 
viously existing  systems  for  such  computations  is'the  capability 
of  including  short-period  information  while  retaining  both  body- 
and  surface-wave  arrivals  on  the  computed  seismograms.  In  this 
context,  by  "short-period"  we  refer  to  the  period  range  covered 
by  the  long-period  instruments  of  the  WWSSN  installations.  In 
practical  terms,  the  successful  accomplishment  of  this  required 
Improvement  is  dependent  upon  highly-optimized  techniques  for 
obtaining  multimode  dispersion-attenuation  information  for  real- 
istic models  of  the  earth,  i.e.,  spherical,  radially  heterogene- 


ous,  anel astic  models. 

Such  optimization  work  has  been  one  of  the  main  interests 
in  our  laboratory  for  several  years.  The  results  of  our  early 
work,  based  on  the  Thomson  (1 950)-Haskel 1 (1953)  technique  and 
on  Knopoff's  (1964)  method  for  treating  flat-layered  structures, 
are  reported  by  Schwab  (1970)  and  Schwab  and  Knopoff  (1970;1971; 
1972;1973).  The  complete  results  of  our  work  on  spheri cal -to- 
f 1 at- s tructure  transformation  techniques,  which  permit  the  use 
of  the  optimized  flat-structure  programs  in  dispersion-attenua- 
tion computations  for  spherical  models  of  the  earth,  are  given  by 
Biswas  and  Knopoff  (1  970),  Schwab  and  Knopoff  (1 971  ;1  972  ;1  973) , 
and  Kausel  and  Schwab  (1973).  In  this  last  reference,  we  have 
also  given  our  initial  outline  of  the  approach  we  have  adopted 
to  handle  the  synthesis  of  multimode  seismograms  once  the  dis- 
persion, attenuation,  source,  and  excitation  functions  have  been 
specified.  This  approach  has  beer,  derived  from  the  basic  theory, 
which  shows  that  the  entire  theoretical  seismogram  for  a dislo- 
cation source  in  a spherical  earth  can  be  expressed  as  a simple 
sum  of  normal  mode  contributions,  as  given  by  Saito  (1967)  and 
Takeuchi  and  Saito  (1972).  Our  first  successes  with  the  system 
we  developed  for  generating  theoretical  seismograms  were  applied 
to  the  interpretation  of  the  seismic  phase  Lg  (Knopoff,  Schwab 
and  Kausel,  1973;  Knopoff  £t  aJL  , 1974).  This  phase  is,  as  we 
demonstrated  in  these  publications,  a multimode  interference 
phenomenon  which  belongs  to  the  surface  wave  portion  of  the 
seismogram.  It  is  only  under  the  present  contract  (see  follow- 

t 

ing  section)  that  we  have  developed  the  system,  which  is  the 


first  we  are  aware  of  for  synthesizing  enti re  seismograms  for 
realistic  models  of  the  earth,  where  both  body  waves  and  sur- 
face waves  appear  on  the  same  record,  and  where  the  period 
content  of  the  record  extends  down  through  the  range  covered 
by  the  long-period  instruments  of  the  WWSSN.  Earlier  work 
of  this  type,  which  was  performed  with  simplified  models  of 
the  earth,  is  summarized  by  Alterman  and  Loewenthal  (1972). 
Sato,  Usami  and  Landisman  (1968)  describe  the  computation  of 
complete  theoretical  seismograms  for  realistic  models  of  the 
earth.  However,  their  results  are  limited  to  ultra  long,  periods 


III.  Objectives,  methods  and  results. 

Data  collection.  To  perform  the  regionalization  of  the  area 
under  investigation  with  the  single-station  method,  we  must  know 
the  orientation  of  the  fault  plane,  depth  of  focus,  and  the  direc- 
tion of  displacement  dislocation  at  the  focus.  These  provide  cal- 
culable corrections  due  to  apparent  initial  phase  of  the  Rayleigh 
wave  signal.  Our  current  Interpretation  procedures,  for  obtaining 
structural  parameters  as  a function  of  depth,  are  based  solely  ci* 
the  phase  velocity  dispersion  of  the  isolated  fundamental  mode; 
hence,  this  correction  can  be  quite  important  in  our  regional iza- 
, tion  studies.  One  problem  which  is  associated  with  accurate  de- 
termination of  the  apparent  initial  phase  is  the  dependence  of 
this  value  on  a precise  knowledge  of  the  structural  parameters 
beneath  the  epicenter,  which  are  used  to  calculate  this  correction. 
To  determine  the  importance  of  this  dependence,  we  have  performed 
a detailed  analysis,  the  results  of  whi*;h  are  reported  by  Frez  and 
Schwab  (1976,  copy  of  page  proofs  appended). 

Our  recent  large-scale  regionalization  work  (Biswas,  1971; 
Kausel  ,1  972;  Leeds, 1973;  Fouda,1  973  ; Kausel  et  a_l_.,1974;  Leeds  e£ 
aj_.  ,1  974;  Knopoff  and  Fouda,1974)  has  shown  the  value  and  ease  of 
| handling  of  seismic  records  of  fundamental  mode  Rayleigh  waves. 

Our  efforts  in  the  present  study  have  therefore  been  focused  on 
this  well-tested  approach. 

The  long-period  records  from  the  47  WWSSN  stations  which  bor- 
der the  region  of  interest  have  been  used  in  the  study.  The  loca- 
tions of. these  stations  are  shown  in  Figure  1.  We  have  also  ob- 
tained several  seismograms  (through  World  Data  Center  B)  of  Soviet 
records  made  in  Central  Asia  on  experimental  long-period  instru- 
ments. However,  the  precision  of  these  recording  plus  an  uncertain 


impulse  response  has  not  permitted  us  to  make  use  of  these  recordings. 

In  order  to  obtain  the  desired  long-period  information  from  WWSSN 
records,  it  is  necessary  to  use  shocks  of  relatively  large  magnitude, 
yet  small  enough  to  be  on-scale  at  the  stations  of  the  network. 
An  example  of  the  intermediate-magnitude  seismicity  of  the  area  is 
given  in  Figure  2.  The  epicenters  are  plotted  for  earthquakes 
having  magnitudes  between  5.9  and  6.6,  which  occurred  during 
the  interval  from  February,  1963  to  February,  1967.  The  choice 
of  this  range  of  magnitudes  is  governed  by  two  considerations. 
First,  experience  has  shown  that  good  long-period  surface  wave 
information  requires  events  with  magnitudes  above  a certain 
value  and,  of  course,  a shock  which  is  so  large  as  to  send  the 
instrument  off  scale  is  useless  for  our  purposes.  Second,  the 
application  of  the  single-station  method  requires  knowledge  of 
the  focal  mechanism.  We  must  therefore  use  events  Terge  enough 
to  allow  us  to  obtain  an  accurate  fault  plane  solution  for  each 
event  we  select  for  processing.  In  addition  to  the  epicenters 
shown  in  Figure  2,  there  are  regions  of  high  seismicity  along 
the  eastern  border  of  the  Kamchatka  peninsula  and  along  the 
Aleutian  arc,  which  may  provide  useful  events  for  this  study. 

With  such  a dense  set  of  stations  around  the  area  to  be 
studied,  and  with  the  regions  of  high,  intermediate  -magnitude  sei smi - 
city  located  within  and  around  the  area,  there  has  been  no 
problem  in  obtaining  sufficient  data  for  the  project.  It  is 
interesting  to  note  that  the  area  is  almost  completely  encir- 
cled: -There  is  only  one  significant  gap  --  between  Japan  and 

Alaska  --  where  stations  do  not  exist.  However,  this  gap  is 

t 

a region  of  high  seismicity  which,  in  the  single-station 


sense*  is  equivalent  to  having  a high  station  density  in  the 
region . 

■ There  are  far  too  many  epicenter-to-statio.n  paths  to 
show  individually,  but  the  limits  of  the  area, which  we  have  cov- 
ered with  a sufficiently  dense  set  of  paths,  are  shown  in  Fig- 
ure 3 by  the  solid  lines.  The  solid  regions  are  those  of 
high,  large-magnitude  seismicity. 

The  first  event  processed  has  the  following  USCGS  speci- 
fication : 

March  7 , 1 966  --  21:29:1  7.4  GMT,  37.3°N,  114. 9°E, 
h = 33  km,  M = 6.0  . 

Our  fault  plane  solution  for  this  event  is  given  in  Figure  4. 

This  first-motion  information  restricts  the  fault  plane  solu- 
tion sufficiently  well  so  that  only  this  information  is  needed 
to  determine  the  required  focal  specification.  The  location 
of  the  epicenter  and  epi center- to-s tation  paths  which  were  se- 
lected for  processing  are  given  in  Figure  5. 

The  computation  of  the  correction  for  the  apparent  ini- 
tial phase  requires  the  specification  of  the  strike  of  the 
fault  plane  4>  , its  dip  6 and  the  direction  of  slip  X , 
and  the  depth  to  the  focus  h.  These  parameters  are  shown  in 

the  fault-plane  geometry  given  in  Figure  6.  Fo  the  above  event, 
our  fault-plane  solution  yields 
4>  = 122°  east  of  north 
' if  82° 

X = 90° 

i 

and  inspection  of  the  Rayleigh-wave  amplitude  spectra  indicates 


( 


a feeal  depth  of 

h = 14.5  km  . 

The  second  event  selected  for  processing  has  the  fol- 
lowing USCGS  specification: 

August  25,  1964  --  13:47:20.6  GMT,  78.2°N,  126. 6°E  , 
h = 50  km,  M = 6.1  . 

The  location  of  the  epicenter,  and  the  epicenter-to-station 
paths  which  were  selected  for  processing  are  given  in  Figure  7. 
The  fault-plane  solution  for  this  event  was  given  by  Sykes 
• (1967)  as 

<J>  = 4°  east  of  north 

’ 6 = 58°  west 

K': 

for  one  possible  fault  plane,  and 

* = 22°  west  of  north 

fe 


! 

f 

6 = 54°  east 

Ik 

if. 

for  the  other. 

Our  result  for  this  event  is 
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4>  = 15°  west  of  north 
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te 
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6 - 58°  west 
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1 A = 260° 

1 

h = 11.5  km  . 

The  angles  <j>  and  6 are  well -constrained  by  the  first-break 
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fault-plane  solution  given  in  Figure  8;  the  additional  informa- 

MU'  , , 

tion  contained  in  the  Rayleigh-wave  amplitude  radiation  patterns 
which  are.  shown  in  Figure  9,  is  required  to  fix  A and  h 

The  epicenter  and  propagation  paths  of  the  third  event 


selected: 


December  26,  1964  --  14:30:29.1  GMT,  51.8°N,  156. 8°E, 
h = 136  km,  M = 5. 7 ( USCGS ) 

are  shown  in  Figure  10,  and  the  results  of  our  first-motion 
fault  plane  analysis  for  this  event  are  given  in  Figure  11. 

The  epicenter  and  propagation  paths  for  the  fourth 
event  selected: 

March  31,  1969  --  07:15:54.4  GMT,  27.7°N,  34.0°E, 
h = 33  km,  M = 6.0  - 6.8 
are  given  in  Figure  12. 

The  epicenter  and  propagation  paths  for  the  fifth 

event  selected:  - 

April  25,  1966  --  23:22:49.3  GMT,  49.3°N,  69.2°E, 

h = 8 km,  M = 5.3 

are  given  in  Figure  20.  The  results  of  the  first-motion  fault 
plane  analysis  for  this  event  are  obtained  from  Zakharova  et. 

al.  (1971): 

<j>  = 55°  west  of  north 
6 = 70° 

1 X - 90° 

The  sixth  event  we  have  used  is  a Chinese  nuclear  explo 
sion  which  has  the  following  USCGS  specification  : 

October  14,  1970  --  07:29:58.6  GMT,  40.9°N,  89.4°E 
The  source  and  propagation  paths  for  this  event  are  shown  in 
Figure  20.  We  assume  the  focal  mechanism  of  the  source  is 
that  of  a simple  explosion  source  at  the  surface  of  the  earth. 

The  combined  Rayleigh  wave  propagation  paths  are  shown 
in  Figure  17,  where  the  continental  areas  are  .emphasi zed,  and 
in  Figure  18,  where  the  Arctic  region  is  emphasized 


Two  events  with  approximately  the  same  epicenter  were 
also  processed.  These  are 

February  6,  1965  - 01:40:33.2  GMT,  53.2°N,  161. 9°W, 
h = 33  km,  M = 6.4  - 6.7 
and 

February  6,  1965  --  16:50:29  GMT,  53.3°N,  161. 8°W, 
h = 33  km,  M = 6.1  - 6.6  . 

The  paths  are  shown  in  Figure  13;  the  fault  plane  solutions 
are  given  in  Figure  14. 

Two  events  with  approximately  the  same  epicenter  were 
the  final  ones  processed.  The  USCGS  specifications  are: 

February  5,  1965  --  09:32:09.3  GMT,  52.3°N,  174. 3°E, 
h = 41  km , M=5.9-6.5 
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February  6,  1965  --  04:02:53  GMT.  52.1“N,  175. 7°E, 
h = 35  km,  M = 5.9  - 6.0  . 

These  paths  are  shown  In  Figure  15;  the  fault  plane  solutions 
are  given  in  Figure  16. 

The  complete  list  of  events  used  in  our  studies  is  as  follows: 
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Two  interesting  points  concerning  the  data-processing 
techniques  have  come  to  light  during  the  present  study.  Both 
involve  the  accuracy  of  the  digitizations  of  the  recorded 
events.  First,  most  of  the  event  records  we  have  used  are 
about  as  large  as  they  could  be  without  going  off  scale.  This 
has  necessitated  a change  in  our  data  processing  techniques 
which  should  be  noted  for  the  information  of  others  involved 
in  this  type  of  work. 

In  the  past,  our  standard  procedure,  when  working  with 
smal 1 er-ampl i tude  recordings,  has  been  to  make  copies  of  events 
from  35  mm  microfilm  records  of  the  WWSSN  seismograms. 

These  copies  are  made  with  a standard  microfilm  reader-printer 
(Itek  18.24  Reader-Printer),  and  the  events  are  then  digitized 
from  these  copies.  We  followed  this  procedure  during  the  ini- 
tial phase  of  the  present  investigation,  but  later  became  con- 
cerned about  the  possibility  of  distortion  in  the  copying  pro- 
cess. Our  tests  comparing  the  phase  velocity  results  obtained 
from  full-size  record  copies  provided  by  NOAA  with  the  results 
obtained  from  our  microfilm  copies  shows  this  concern  to  be 
valid.  Our  conclusion  is  that,  when  working  with  large-ampli- 
tude recordings  sich  as  those  employed  in  the  present  study, 
one  must  use  full-size  record  copies;  one  must  not  use  a micro- 
film copier  which  forces  the  analysis  to  use  several  prints, 
spliced  together,  to  form  the  record  copy  from  which  the  digi- 
tization is  obtained. 
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The  specific  period  ranges,  which  were  used  in  the  present  in- 
version, are  shown  in  Table  1. 

The  phase  velocities  from  the  five  earthquakes  and  one 
nuclear  explosion  for  the  32  paths  in  question  (Figure  20)  ap- 
pear to  assort  themselves  into  two  groups  (Figures  21  and  22). 

The  paths  with  higher  phase1  vel oci ties  are  generally  those 
that  cross  the  stable  platforms  and  shields  (such  as  paths  from 
the  Hsingtai  earthquake  to  Scandinavian  (KEV)  and  German  (STU) 
stations;  typical  of  the  lower-velocity  group  is  the  phase  velo- 
city on  the  paths  from  the  Red  Sea  earthquake  to  the  southern 
Asiatic  stations  (SHL,  MAN,  etc.)  . Phase  velocities  for  typi- 
cal shield  regions  and  young  stable  continental  regions  (Biswas 
and  Knopof f , 1 974 ) are  shown  for  comparison.  The  incompatibility 
of  most  of  the  phase  velocity  observations  for  Eurasian  paths 
with  the  "standard"  curves  for  homogeneous  regions  testify  to 
the  inhomogeneous  nature  of  the  Eurasian  continent. 

The  first  step  in  the  structural  analysis  required  the  divi- 
sion of  the  area  into  subregions,  in  each  of  which  the  structural 
'parameters  are  assumed  to  be  laterally  homogeneous.  Our  first 
subdivision  is  shown  in  Figure  23,  and  was  based  on  the  use  of 
the  worldwide  tectonic  map  (Khain  and  Muratov , 1 969) . Parts  of 
this  map  are  too  detailed  for  our  purposes;  however,  they  provide 
the  basis  for  our  regionalization.  We  have  subdivided  the  Eur- 
asian area  into  six  regions  which  represent  structural  averages 
of  a broad  nature.  In  the  subdivision,  we  have  been  guided  by 
the  development  of  continental  regionalization  based  on  inversion 
of  surface  wave  phase  velocity  data  (Knopoff ,1 972) . We  have  used 
the  tectonic  and  heat-flow  maps  to  indicate  where  certain  typical 


geophysical  provinces,  as  expected  from  surface  wave  studies, 
are  to  be  found.  The  regions  we  have  used  are 

1.  Ancient  PreCambrian  Shields 

2.  Younger  (seismical ly)  Stable  Platforms 

3.  The  Himalayan-Alpide  Mountain  Belt 

4.  The  Tibetan  Plateau 

5.  The  Si nki ang-Mongol i an  Seismic  Zone 

6.  The  Chinese  "Stable"  Region 

The  ancient  PreCambrian  Shields  are  well-defined,  both  as. 
tectonic  (from  the  map)  and  seismic  units.  The  same  comments 
apply  to  the  Younger  Stable  Regions.  Although  few  seismic  data 
have  been  accumulated  to  date  with  regard  to  folded  mountain 
belts  --  the  exception  is  the  Alps  (Knopoff  et  al  ,1966)  - for  convenience 
in  computation,  we  establish  this  zone  as  a single  seismic 
province.  Region  4,  the  Tibetan  Plateau,  is  definable  in  terms 
of  the  enormous  elevation  associated  with  this  region.  We  have 
no  justification  for  setting  up  regions  5 and  6,  except  that 
these  seem  to  be  regions  definable  in  terms  of  seismic  activity, 
the  modern  seismic  activity  of  these  regions  is  inconsistent 
with  surface  geology  (very  old  rocks,  in  some  cases)  and  to  us 
indicates  the  presence  of  some  anomalous  behavior  in  the  upper 

• 

mantle.  Although  these  regions,  within  themselves,  encompass 
widely  differing  geologic  structures  and  widely  varying  seismic 
activity,  we  have  assumed  that  each  of  these  regions  is  homo- 
geneous. In  the  inversion,  we  have  not  assumed  a priori,  that 
these  two  regions  have  a structure  which  is  similar  to  that  of 
any  of  the  other  numbered,  regions;  should  the  inversion  give  the 
result  that  regions  5 and  6 are  the  same  as  one  of  the  other  re- 
gions, this  can  be  used  to  simplify  the  inversion  in  a later 
stage.  It  can  be  seen  from  th‘>  table  of  period  ranges  of  the 
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observations  that  the  inability  to  obtain  long-period  informa- 
tion,from  the  observations  of  the  Tashkent  earthquake  and  the 
Lop  Nor  nuclear  explosion,  limits  our  ability  to  resolve  deeper 
structure  under  region  4.  This  will  be  more  evident  in  the  re- 
sults of  the  inversion  to  be  described  below. 

Two  small  regions  are  not  included  in  our  inversion:  The 

South  China  Sea,  and  the  Sea  of  Okhotsk.  For  the  South  China 
Sea,  we  have  assumed  the  structure  to  be  known,  and  to  be  that 
for  typical  marginal  seas.  The  phase  delays  for  this  region 
• are  taken  from  two  single-station  phase  velocity  observations 
in  marginal  seas  obtained  by  Leeds  (1973)  and  from  two  phase 
velocity  determinations  by  the  two-station  method  obtained  by 
Knopoff  (unpublished)  across  the  Philippine  Sea.  Phase  velocity 
. corrections  for  the  Sea  of  Okhotsk  were  obtained  theoretically 
from  a structure  derived  from  Kosminskaya  et  al . (1969),  in 
which  the  Okhotsk  depression  has  a 25  km  crustal  structure; 
we  have  used  an  oceanic  mantle  below  this  crust.  These  phase 
corrections  have  been  applied  to  travel  times  for  those  paths 

< 

that  traverse  these  two  regions. 

The  total  path  length  summed  over  all  event-paths  in  each 
region  is  given  in  the  following  table; 


Region 

Path 

Length  (km) 

TABLE  2 

Percentage  of  Total 
Path  Length 

Weighted  Percentage  of 
Total  Path  Length 

1 

20970 

12.1 

13.0 

2 

.'53150 

30.6 

30.4 

3 

38411 

22.1 

22.5 

4 

15292 

8.8 

i 

6.9 

5 

35994 

20,7 

20.7 

6 , 

9979  , 

5.7 

6.5 

TOTAL 

173796 

100.0 

100.0 

Since  the  harmonic  analysis  for  phase  has  not  been  carried  out 
over  the  same  band  of  periods,  a column  of  the  above  table  has 
been  included  to  give  the  product  of  path  length  by  number  of 
period  estimates  of  phase.  By  either  method  of  estimation,  the 
very  low  fraction  of  sampling  in  regions  4 and  6 will  give  us 
large  uncertainties  in  the  structure  obtained  from  the  inversion 
for  these  regions. 

We  have  inverted  the  data  under  the  assumption  that  a sim- 
ple ray  theory  for  surface  waves  applies.  That  is,  the  phase 
shift  for  a seismic  wave  passing  through  a given  region  is  com- 
puted as  though  that  region  were  laterally  infinite  in  extent 
and  uninfluenced  by  the  presence  of  neighboring  regions,  no  mat- 
ter how  close  the  great  circle  path  approaches  the  regional  boun- 
daries. This  assumption  is  evidently  untenable  but  provides  a 
basis  for  starting  an  inversion;  future  work  is  planned  to  test 
the  accuracy  of  ■'his  assumption  and  to’  provide  a means  of  approx- 
ima te  correcti on . 

The  inversion  proceeds  as  in  the  method  described  by  Leeds 
et  al  . (1  974).  We  calculate  the  phase  travel  time  for  the  ntfl 
path  and  the  pth  period  as 

6 

= y z.  s . 

nP  i=i  in  ip 

where  Z]r]  is  the  path  length  of  the  nth  path  in  the  1th 
region  (i=l,...,6)  and  s^  is  the  (theoretical)  phase  slow- 
ness for  the  ith  region  at  the  pth  period.  The  phase  slow- 
nesses s.jp  are  functions  of  the  model  parameters  in  each  region. 


The  more  model  parameters  we  take  in  each  region,  the  greater 
the  uncertainties  in  the  determination  of  the  parameters.  We 
minimize 


l £ ( t - tnn)2 
n p nPo  rip' 


where  tnp0  is  the  observed  travel-time  for  the  nth  path  at 
the  pth  period.  The  minimization  takes  place  with  rdspect  to 
the  choice  of  model  parameters. 

, The  inversion  becomes  less  and  less  accurate,  i.e.,  we  ob- 

tain large  variances  in  the  model  parameters,  as  the  number  of 
model  parameters  increases.  It  would  be  pleasing  to  be  able  to 
solve  for  the  properties  of  the  crust  in  each  of  the  six  regions. 
However,  this  a)  requires  much  precise  data  at  periods  shorter 
than  30  sec.,  b)  increases  the  number  of  model  parameters  sig- 
nificantly, and  c)  undoubtedly  pushes  our  postulate  of  lateral 
homogeneity  in  each  of  the  regions  to  an  untenable  extreme.  We 
have  therefore  used  model  crusts  for  each  of  the  regions  which 
a')  seem  to  be  plausible  when  compared  with  results  for  other  sim- 
ilar parts  of  the  earth  where  ob:  o vations  exist  (such  as  loca- 
ations  typical  of  regions  1,  2 and  3),  and  b)  agree  with  Soviet 
refraction  data  when  available.  When  large  residuals  were  en- 
countered at  short  periods,  such  as  in  the  case  of  region  3 ( A 1 - 
p i de-Hima 1 ayan  belt)  and  region  4 (Tibetan  plateau),  we  were  ob- 
liged to  Introduce  more  low-velocity  material  into  the  crust. 

This  was  done  by  keeping  crustal  velocities  fixed  and  increasing 
crustal  thickness.  In  these  two  cases,  this  procedure  leads  to 
extraordinarily  thick  crusts.  It  should  be  realized  that  these 
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model  crustal  thicknesses  are  consequences  of  the  procedures 
used;  if  we  had  chosen  to  lower  the  crustal  velocities,  the 
thicknesses  would  have  been  less. 

We  have -pa ranieter i zed  the  mantle  into  a lid,  channel  and 
subchannel;  each  of  these  layers  in  a given  region  is  homogen- 
eous. The  subchannel  region  terminates  at  a depth  of  420  km. 
Below  this  depth,  we  place  a standard  lower  mantle  platform 
under  all  regions.  Once  again,  in  order  to  reduce  the  number 
of  degrees  of  freedom  in  the  inversion,  we  have  fixed  the  lid 
S-wave  velocity  at  4.65  km/sec  and  the  subchannel  S-wave  velocity 
at  4.8  km/sec.  The  value  of  4.65  km/sec  for  the  lid  was  chosen 
because  it  arises  frequently  in  inversions  for  many  other  parts 
of  the  world.  In  the  most  obvious  case  in  which  a 4.65  km/sec 
lid  was  not  observed,  namely  for  the  Western  United  States  (Bis- 
was and  Knopoff ,1 974)  in  which  the  subcrustal  material  has  S- 
wave  velocity  around  4.3  km/sec,  we  can  assume  that  a model  with 
a 4.65  km/sec  lid  with  zero  lid  thickness  was  accepted  by  the 
inversion  and  that  the  4.3  km/sec  value  is  representative  of 
channel  material  rising  almost  to  the  crust.  Should  the  lid 
I velocity  in  some  region  be  less  than  4.65  km/sec  in  reality, 
then  crustal  thicknesses  can  be  reduced;  we  suspect  that  for 
regions  3 and  4 this  possibility  may  exist.  A similar  comment 
can  be  made  about  the  subchannel  velocity:  if  the  subchannel 

velocity  should  turn  out  to  be  less  than  4.8  km/sec  in  some  re- 
gions, channel  thicknesses  will  be  reduced. 

The  parameterization  thus  includes  only  two  adjustable 
model  parameters  for  each  region.  These  are  the  lid  thickness 
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and  the  channel  S-wave  velocity.  A thirteenth  parameter  is 
the  subchannel  thickness,  which  is  presumed  to  be  uniform  across 
the  entire  continent  and  hence  has  the  same  value  under  each 
region.  Since  the  crustal  thickness  and  the  depth  to  the  420  km 
interface  are  fixed,  the  parameterization  of  lid  and  subchannel 
.thicknesses  is  equivalent  to  a parameterization  of  the  depth 
below  the  surface  of  the  top  and  bottom  of  the  channel.  This 


. parameterization  has  13  degrees  of  freedom. 

After  adjustment  of  the  crust  by  the  procedures  described 

above  (with  interpretation  of  crustal  parameters  according  to 
the  caveats  expressed  above),  the  parameterization  and  cross- 
sections  used  in  a linearized  inversion  procedure  are  given  in 
Table  3.  The  superficial  sedimentary  layer  that  is  introduced 
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TABLE  3 

Model  parameters  for  inversion  of  phase  velocity  data 

Thickness  Depth  - f3(km/sec)  ot(k(n/sec)  p(gm/cm3) 

— 0 

CRUST  . (different  crustal  models  for  each  region) 

(fixed) 

LID  VAR  4.65  8.17  3.45 


CHANNEL  VAR  7.80  3.45 


SUBCHANNEL  VAR  4.80  8.80  3.65 


420( fixed) 


TABLE  3 (cont'd) 

Detailed  Crustal  Models  Used  in  Inversion 
(thickness  in  km,  velocities  in  km/sec,  density  in  gm/cm 
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The  starting  set  of  thirteen  parameter  in  the  linearized 
inverse  is  (in  the  usual  units): 


TABLE  4 


Region 

6CH 

hLID 

hSUB 

1 

4.51 

no 

2 

4.39 

113 

3 

4.30 

54 

If 

>0 

4 

4.29 

90 

5 

4.08 

92 

6 

4.38 

103 

' 

Depth  to 
top  of 
channel 

155 

158 

119 

167 

137 

148 


Depth  to 
bottom  of 
channel 


270 


The  thirteen  parameters  in  the  rectangular  box  in  Table  4 are 
those  varied  in  the  inversion. 

The  accuracies  of  the  measurements  were  taken  to  be  the 
same  as  those  used  by  Leeds  et  al.  (1974),  namely  a measurement 
error  presumed  to  be  7 sec  for  periods  less  than  70  sec  and  one- 
tenth  the  period  for  periods  greater  than  70  sec  (hence,  at  a 
period  of  357  sec,  the  presumed  standard  error  is  36  sec).  • 

In  the  inversion,  an  iteration  process  has  been  used  in 
which  the  matrix  of  partial  derivatives  6 was  recalculated 
whenever  we  moved  into  a new  portion  of  parameter  space. 

This  process  was  continued  until  we  obtained  satisfactory  con- 
vergence of  the  variable  model  parameters.  The  thirteen  eigen- 
values to  the  product  matrix  of  partial  derivatives  GTG  , which 
were  obtained  in  the  final  stage  of  the  Iteration  process,  are 
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given 

in  Table  5.  Each  eigenvalue  corresponds  to  an  eigenvector 

■\o.^V 

uv . 

which 

in  every  case 

points  in  a direction  c‘ose  to  one  of  the 

: Jpj1  •■ 

;fe- 

§■ 

;$>■■ 

mVJv' 

thirteen  parametric 

axes.  Thus,  each  eigenvector  can  be  said 

to  be 

a discriminant 

for  one  of  the  thirteen  degrees  of  freedom 

1^ 

in  the 

model . From 

the  theory  of  the  inversion  analysis,  we 

fe: 

tfe: 

know  that  the  eigenvectors  corresponding  to  the  largest  eigen- 
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give  reliable 

structural  information  concerning  the  model 
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its  corresponding  model  parameter. 
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TABLE  5 

1 

1 

1 '■ 

Eigenvalue 

Model  parameter  most  closely  resolved 
by  eigenvector 

1. 

13.68 
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2. 

8.43 
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From  this  table,  it  can  be  seen  that  in  general  the  data  con- 
tain the  greatest  amount  of  information  about  the  channel  velo- 
cities and  the  least  information  about  the  lid  thicknesses.  In- 
general, the  data  contain  more  information  about  regions  2,  3 
and  5 than  the  other  three  regions.  This  is  due  to  the  higher 
percentage  of  total  path  length  which  samples  regions  2,  3 and 
5.  (For  complete  details,  see  Table  2 and  Fig.  24). 

The  variance  of  each  model  parameter  decreases  as  succes- 
sive models  approach  the  final  one,  which  demonstrates  the  con- 
vergence of  our  inversion  procedure.  Fig.  25  gives  our  final 
results:  the  best  model  (in  a "least-squares"  sense)  and  the 

corresponding  standard  deviations  of  each  model  parameter. 

The  quality  of  the  agreement  of  our  model  to  the  experimental 
data  is  expressed  by  the  parameter  R : 

[\  n y2 


where  n is  the  total  number  of  observations,  A t ^ is  the 
travel  time  residual  and  o.  is  the  presumed  standard  error 
in  the  data.  If  the  error  for  each  observation  has  been  cor- 
rectly estimated,  R should  be  about  1.0  for  the  final  model. 
A value  of  R<1  suggests  that  the  errors  are  smaller  than  es- 
timated; R>>1  indicates  a poorly-fitting  model.  The  value 
of  R is  1.126  for  the  final  model  from  our  inversion. 


From  Fig.  25,  we  also  see  that  the  upper  mantle  structure 
for  regions  I and  II  are  very  similar.  The  lid  in  region  I is 
somewhat  thicker.  The  channel  shear-wave  velocity  for  region  I 
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is  slightly  less  than  that  in  region  II,  but  the  uncertainty 
of  the  model  parameter  in  region  I is  rather  large;  thus,  the 
existence  of  a low-velocity  channel  in  region  I is  uncertain. 
The  most  striking  result  of  this  inversion  is  the  very  thin 
lid,  and  moderate  shear-wave  velocity,  in  the  channel  for  re- 
gion III  and  V.  These  two  regions  exhibit  tectonic  activity 
and  high  seismicity.  Region  IV  is  that  in  which  the  Eurasian 
and  Indian  continents  are  in  collision.  Although  the  upper 
mantle  structure  for  this  region  appears  to  be  very  similar  to 
those  of  regions  I and  II,  the  uncertainties  in  the  model  para- 
meters are  large.  This  is  due  1)  to  the  low  percentage  of 
total  path  length  in  this  region,  and  2)  to  the  fact  that  most 
of  the  paths  which  sample  this  region  have  only  short-period 
information.  Region  VI  has  a thick  lid  and  a pronounced  low- 
velocity  channel,  but  the  standard  deviations  are  rather  large; 
additional  paths  will  be  required  to  improve  the  resolution  of 
the  structural  parameters  in  this  region. 

A number  of  control  experiments  and  additional  measure- 
ments are  planned  for  our  next  contract.  These  include: 

1)  Acquisition  of  additional,  longer-period  phase  velocity  in- 
formation which  samples  the  regions  not  wel 1 -resol ved  in 
our  current  inversion.  These  regions  include  Tibet  and 
Southeastern  China. 

2)  Investigating  the  stability  of  the  results  by  varying  the 
location  of  the  provincial  boundaries.  Some  of  this  work 
has  already  been  done,  but  more  needs  to  be  done  in  this 


regard.  • 

3)  Investigating  the  corrections  to  the  raw  data  to  be  applied 





z 


to  take  into  account  events  whose  geometrical  paths  are 
close  to  and  parallel  to  provincial  boundaries. 

4)  Seeing  whether  it  is  possible  to  reduce  the  number  of 
degrees  of  freedom  still  further  beyond  those  already 
considered,  possibly  by  reducing  the  number  of  provinces, 
and  by  constraining  the  bottom  of  the  lid  to  fit  the  melt- 
ing curve  for  peridotite. 

5)  To  acquire  Love  wave  data  plus  SRO  and  high-gain  data  for 
this  region.  The  SRO  data  will  be  especially  valuable  in 
helping  us  to  resolve  the  deeper  structure  in  these  regions. 

6)  Investigating  the  linearity  of  the  solution  by  varying  the 
starting  model . 


36. 


Arctic  Study.  We  have  measured  fundamental  mode  Rayleigh 
waves  over  a number  of  paths  crossing  the  Arctic  Ocean.  We 
have  used  as  sources  four  earthquakes  whose  focal  parameters 
are: 


Lena  River 

Aug. 

25, 

1964 

13:47:20.6 

78 . 2 °N 

126. 6°E 

Eastern  Aleutians(l) 

Feb. 

6, 

1965 

01  : 40 : 33 . 2 

53 . 2°N 

161 .9°W 

Eastern  A1 eutians (2) 

Feb. 

6, 

1 965 

16:50:23.6 

53. 3°N 

161 . 8°W 

Western  Aleutians 

Feb. 

6, 

1965 

04:02:53 

52.1 °N 

175. 7°E 

For  each  of  these  earthquakes  we  have  obtained  initial  phases 
either  from  the  fault  plane  solution  or  from  the  radiation  pat- 
tern for  Rayleigh  waves  (as  in  the  case  of  the  Lena  River  dis- 
cussed above).  We  have  obtained  phase  velocities  by  the  single- 
station method  for  seven  paths  crossing  the  Arctic  over  the 
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It  can  be  seen  that  none  of  these  are  purely  oceanic  paths. 
The  heavy  line  outlines  our  estimate  of  the  boundary  between 
the  continental  shelf  and  the  deep  ocean  basins.  The  fraction 
of  the  geometrical  path  that  each  event  has  in  the  oceanic  part 
is  as  follows: 


Event 


Total  path  length  Oceanic  length  Fraction  oceanic 


1 . Lena  - ESK 

4816  km 

1904  km 

.40 

2.  Lena  - KTG 

3386 

1823 

.54 

3.  East  Aleut(l)  - K0N 

7483 

4603 

.62 

4 . East  Aleut(l ) - KTG 

. 5933 

916 

.15 

5.  East  A1 eut (2 ) - KEV 

6347 

2334 

.37 

f).  East  Aleut(2)  - ESK 

7818 

i 

2714 

.35 

7 . West  Aleut.  - KEV 

6263 

844 

.13 

, ...  '■  ^ 
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To  reduce  the  available  data  to  information  regarding 
purely  oceanic  events,  we  have  decided  to  use  the  phase  veloci- 
ty data  for  the  Lena  River  event  recorded  at  KEV  (see  discus- 
sion above  for  Eurasia)  as  a typical  continental  value  and  to 
subtract  these  numbers,  for  the  appropriate  path  length  con- 
tribution, from  the  phase  delays  observed  for  the  above  7 
path-events.  Unfortunately,  the  two  events  East  Aleut. (1)  - KTG 
and  West  Aleut. -KEV  have  such  small  parts  of  this  total  path 
that  are  oceanic  that  we  are  subtracting  two  numbers  of  com- 
parable size  and  the  result  is  quite  unstable.  The  unreliabil- 
ity of  the  reduction  results  for  these  two  cases  has  obliged 
us  to  exclude  them  from  our  data  set.  Accordingly,  we  have  in- 
vestigated the  implications  of  the  phase  velocity  results  for 
the  five  remaining  paths.  The  relevant  data  are  given  in  the 
following  Table: 


(Pure)  Oceanic  Phase  Velocities  (km/sec) 


T(sec) 

Lena-ESK 

Lena-KTG 

E.A1 (1 ) - K0  N 

E . A1 ( 2 ) -KEV  E. 

. A1 ( 2) -ESK 

208 

4.64 

(4.83) 

4.62 

4.55 

4.52 

192 

4.52 

(4.61) 

4.51 

4.42 

4.42 

167 

4.38 

4.32 

4.38 

4.22 

4.28 

139 

4.21 

4.14 

4.24 

4.09 

4.15 

119 

4.08 

4.06 

4.12 

4.04 

4.11 

100 

4.01 

3.99 

4.08 

3.96 

4.07 

69. 

3.95 

‘3.90 

4.00 

3.91 

4.00 

50 

3.91 

3.82 

3.98 

3.89 

3.92 
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the  deep  Arctic.  The  best  we  can  do  Is  consider  the  deep  Arctic 
as  a single  province  and  investigate  the  consequences  of  inver- 
ting an  "average"  phase  velocity  for  the  region.  The  average 
phase  velocity  is  obtained  from  the  above  table  by  weighting  by 
the  oceanic  path  length  In  each  case.  The  result  Is  (omitting 
the  quantities  in  parentheses): 


T(sec) 

c(km/sec) 

*$?*.- 

208 

4.59 

192 

4.47 

167 

4.32 

■ 

139 

4.18 

ft’ 

m- 

V: 

119 

4.10 

rk' 

1 

100 

4.03 

f 

69  - 

3.97 

1; 

50 

3.92 

These  results  can  be  compared  with  those  obtained  for  Pacific 
paths  by  Leeds  (1973)  from  inversion  of  trans-Pacific  phase 
velocity  data  by  methods  similar  to  those  described  above  for 
trans-Eurasian  paths.  After  inversion  and  determination  of  the 
, pure-age  cross-sections,  the  pure-age  phase  velocities  for  the 
Pacific  and  be  derived;  these  are  shown  in  Fig. 27  for  Pacific 
ages  0-10  my,  20-40  my,  85-110my.  The  Arctic  data  points  are 
shown  as  circles.  The  Arctic  cross-section  averages 

out  to  about  a 30  my  Pacific  structure.  According  to  the  model 
of  Parker  and  Oldenburg  (1973),  the  lithospheric  thickness  Is  a 
function  of  the  age  z(t)  . 9.4t*  km  .with  t the  age  in  my. 
Thus,  assuming  that  the  lid  and  channel  velocities  are  of  the  same 
order  as  In  the  Pacific,  the  age  at  which  the  Arctic  began  to  open 
Is  calculated  to  be  about  70  my  (before  present),  an  unexpectedly 

small  quantity. 
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Theoretical  seismograms.  Prior  to  the  availability  of 
the  final  regionalization  from  the  first  part  of  our  investi- 
gation, the  thrust  of  our  v/ork  with  time  series  synthesis 
has  been  directed  toward  improving  the  efficiency  of  existing 
computational  techniques.  This  improvement  is  required  to 
permit  us  to  extend  the  information  contained  on  the  theoreti- 
cal seismograms  down  through  the  period  range  covered  by  the 
long-period  instruments  of  the  WWSSN.  Although  up  to  the  pre- 
sent time  we  have  concentrated  on  laterally  homogeneous  struc- 
tures, in  all  other  respects  our  models  of  the  earth  have  been 
highly  realistic:  approximately  200  layers  were  used  to  model 
the  radial  heterogeneity  of  the  crust-mantle  system  of  a spher- 
ical earth,  and  the  intrinsic  anelasticity  was  handled  in  an 
exact  manner. 

A summary  of  the  general  methods  we  have  applied  in  our 
computations  is  given  by  Kausel  and  Schwab- (1973)  and  Knopoff 
crt  aJL  (1  973).  An  elaboration  of,  and  certain  justifications 
for,  these  procedures  have  recently  been  given  by  Schwab  and 
Kausel  (1  976  ).  A recent  contribution  by  Calcagnile  ejt  al_.  ( 1 976  , 
preprint  appended),  completed  under  this  contract,  is  also  pertinent 
here,  when  only  the  surface-wave  portion  of  the  theoretical 
seismogram  is  desired. 

Algorithmic  development  and  programming  improvements, 
which  were  carried  out  under  the  present  contract,  are  contained  in 
several  appended  reprints  and  preprints:  Schwab  et^  al_.  (1974), 

Nakanishi  et  al_.(1  975),  Nakanishi  et  al_.  (1  976a),  and  Nakanishi 
et  a!..  (1  976b).  ; 


new  work  is  being,  and  will  to  be,  carried  out  under 
our  next  contract,  and  involves  our  preliminary  attempts  to  com- 
pare the  computed  seismograms  with  the  entire  record  from  long- 
Period  instruments  of  the  WWSSN,  The  initial  results  of  this 
are  reported  by  Kausel,  Schwab  and  Mantovani  (1976,  preprint 
appended),  and  by  Mantovani  et  aj,  (1976a, preprint  appended).  0, 

further  manuscript  along  these  lines  (Mantovani  et_ji,  ,976b)  is 
approaching  completion. 


Our  further  plans  with  this  work  involve  the  addition 
of  more  of  the  higher  modes  in  the  seismogram  synthesis.  Cur- 
rently, we  are  using  11  modes;  we  would  like  to  increase  this 
to  21.  Also,  we  are  now  ready  to  make  the  first  approach  to 
the  treatment  of  lateral  heterogeneity.  This  will  involve  the 
use  of  three  separate  structures;  one  to  model  the  parameters 
beneath  the  epicenter,  and  to  be  used  in  the  generation  of  ini- 
tial excitation  and  phases  at  the  source;  one  structure  as  an 
"average"  over  the  prop  'ation  path;  and  a final  structure  to 


model  that  beneath  the  instrument.  Finally,  to  date  we  have 
limited  our  work  with  theoretical  seismograms  to  the  SH  compo 


nent  of  motion;  we  are  now  interested  in  extending  the  treat- 
ment to  include  the  P-SV  (Rayleigh)  components. 

The  final,  practical  purpose  of  our  work  under  this  con- 
tract --  application  of  our  results  to  the  discrimination  prob- 
lem - involves  comparison  of  theoretical  and  experimental 
seismograms.  It  is  therefore  Important  that  we  have  as  accurate 
a means  as  possible  of  obtaining  the  instrumental  constants  from 
the  impulse  response  of  the  experimental  record.  These  con- 
stants then  permit  us  to  include,  with  a minimum  of  error,'  the 


instrumental  response  on  the  theoretical  seismogram.  Our  im- 
proved scheme  for  inversion  of  the  impulse  response  to  obtain 
the  electromagnetic  instrument  constants,  which  has  been  de- 
vised under  the  present  contract,  is  described  by  Mitchel, 

Knopoff  and  Schwab  (1975,  reprint  appended)  and  by  Mitchel  (1976, 
preprint  appended ) . 
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V.  Figure  Captions 


m.  I.  Locations  of  HHSSN  stations  and  new,  high-gain  instill.- 
tions  to  be  used  in  the  investigation. 


FIG.  2. 


FIG.  3. 


FIG.  5. 


FIG.  G . 


Example  of  tfio  large-magnitude  seismicity  of  the 
region  we  propose  to  study.  The  regions  of  high, 
large-magnitude  seismicity  along  the  eastern  border 
of  the  Kamchatka  peninsula  and  along  the  Aleutian 
arc  may  also  provide  useful  events  for  the  study. 


Limits  (solid  lines)  of  the  region  we  propose  to  cover 
with  a dense  set  of  epicenter-to-station  paths.  Solid 


- S ens  are  the ■' "■  Ma1' 
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Sc  magnitude  seismicity. 


FIG.  4.  First-motion  information  and  fault  plane  solution 

for  event  occurring  at  21:29:17.4  GMT  on  March  7, 
1966. 


Location  of  epicenter  and  epicenter-to-station 
paths  selected  for  processing  from- event  occurring  at 
21 : 29 : 1 7 . 4 GMT  on  March  7,  1966. 


Geometry  an.d  coordinate  systems  of  fault  plane, 
focus,  and  epicenter. 
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FIG.  7.  Location  of  epicenter  and  cpicenter-to-statlon 

paths  selected  for  processing  for  event  occurring 
at  13:47:20.6  GMT  on  August  25,  1964. 

FIG.  8.  First-motion  Information  to  be '.used  for  fault-plane 
solution  for  event  occurring  at  13:47:20.6  GMT  on  August  25, 
1964.  The  solution  given  by  Sykes  (1967)  is  Indicated  by  dotted 

FIG.  9.  Rayleigh-wave  amplitude  information  to  be  used  for 

fault-plane  solution  for  event  occurring  at  13:47:20. 
GMT  on  August, 25,  1964.  The  central  set  of  radiation 
patterns  are  the  results  of  theoretical  computations 
based  on  the  fault  plane  solution  given  In  the  text. 

The  other  four  radiation  patterns  depict  the  expert- 
men ta 1 res  u 1 ts  . 

FIG.  TO.  Location  of  epicenter  and  epicenter-to-statlon  paths 


.FIG. 


FIG. 


FIG. 


selected  for  processing  from  event  occurring  at 
14:30:29.1  GMT  on  December  26,  1964. 

n.  Results  from  first-motion  fault  plane  analysis  for 

event  occurring  at  14:30:29.1  GMT  on  December  26,  1964 

12.  Location  of  epicenter  and  epicenter-to-statlon  paths 
selected  for  processing  for  event  occurring  at 
07:15:54,4  GMT  on  March  31,  1969, 

13.  Location  of  epicenters  and  epicenter-to-statlon  paths 
• selected  for  processing  for  events  occurring  at 

01:40:33.2  (dashed  lines)  and  16:50:29  GMT  (solid 
lines)  on  February  6,  1965. 
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FIG.  14.  Results  from  first-motion  fault  plane  analysis 

for  events  occurring  at  01:40:33.2  (Fig. 14a)  and 
16:50:29  GMT  (Fig. 14b)  on  February  6,  1965. 

FIG.  15.  Location  of  epicenters  and  epicenter-to-station 
paths  for  events  occurring  at  09:32:09.3  GMT  on 
February  5,  1965  (dashed  line)  and  at  04:02:53 
GMT  on  February  6,  1965  (solid  line). 

FIG.  16.  Results  from  first-motion  fault  plane  analysis 
for  events  occurring  at  09:32:09.3  GMT  on  Febru- 
ary 5,  1965  and  at  04:02:53  GMT  on  February  6, 

1 965. 

FIG.  17.  Combined  Rayleigh  wave  propagation  paths.  Con- 
tinental results  are  emphasized  here. 

fiG.  i o . Cuiiiuined  Rayieigii  wave  propagation  paths.  Arctic 
results  are  emphasized  here. 

FIG.  19.  Sample  phase  velocity  results  for  the  paths  noted. 
The  paths  are  shown  in  Figures  5 and  7. 
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